The chiral condensate in holographic models of QCD 
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Bottom-up holographic models of QCD, inspired by the AdS/CFT correspondence, have shown 
a remarkable degree of phenomenological success. However, they rely on a number of bold as- 
sumptions. We investigate the reliability of one of the key assumptions, which involves matching 
the parameters of these models to QCD at high 4D momentum and renormalization scale fi^. 
We show that this leads to phenomenological and theoretical inconsistencies for scale-dependent 
quantities such as (qq). 
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There are still no systematic analytic tools to study the 
strong-coupling dynamics of QCD, except for models that 
probe certain limited classes of observables. Thus, for 
example, one can use chiral perturbation theory for some 
low-energy observables, but for more general ones one is 
forced to resort to more phenomenological approaches 
such as chiral soliton models. In recent years, 'bottom- 
up' holographic models of QCD have emerged as another 
approach to the low energy phenomen olog y of QCD, and 
have attracted considerable interest [l|-[l9{ . 

In these models, QCD in the large Nc limit is taken to 
be dual to a classical 5D theory in a curved space, and the 
parameters of the 5D model are matched to their corre- 
sponding values in large Nc QCD[23,[1H, with the field 
content of the 5D models chosen to match the low en- 
ergy chiral symmetry of QCD. In contrast to approaches 
like chiral perturbation theory, these models allow the 
computation of meson spectra and couplings, at least in 
principle. Even very simple 5D models seem to show a 
remarkable phenomenological success when compared to 
data (see, for instance, ref. [EB)- 

There are, however, some phenomenological puzzles 
with the models. For example, quark masses as extracted 
by ref. [l| and ref. Q are different by a factor of four. 
This inconsistency, which we address in this Letter, is in 
fact symptomatic of issues which are more serious than 
they may first appear. The resolution of these issues, as- 
sociated with the treatment of operators which are scale 
dependent in QCD in holographic models, will shed light 
on the the matching of the 5D models to QCD. 

Holographic models of QCD are motivated by the con- 
jectured dualities between some types of gauge theories 
and higher-dimensional gravity theories. The most well- 
known such duality is the AdS/CFT correspondence [H- 
[23 |. where a conformal field theory, J\f = A SU{Nc) su- 
per Yang-Mills theory in the large Nc limit, is dual to a 
type IIB string theory on AdS^ x S^, where AdSs is 5D 
Anti-DeSitter space, and is the 5-sphere. When the 't 
Hooft coupling of the field theory is large, the AdS^ x 
physics is described by weakly-curved classical supergrav- 



ity. The CFT lives on the boundary of AdS^ , in the sense 
that each operator O of the CFT is identified with a bulk 
field in AdS^ x S^, the boundary values of the bulk 
fields (po' are the sources for the CFT operators, and the 
supergravity partition function Zsq is identified as the 
generating functional of the CFT correlation functions: 
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CFT correlation functions in the strongly-coupled do- 
main can be calculated by evaluating the supergravity 
action Ssg on the classical solution for fields (fP that 
approaches a specified boundary value ( t>n-, and taking 
functional derivatives with respect to (j)^ [23, [23| ■ 

The bottom-up approach to holographic QCD gener- 
ally consists of guessing a 5D background and field con- 
tent that captures some aspects of large Nc QCD for 
some observables of interest. In AdS/CFT, the confor- 
mal invariance of the CFT corresponds to the fact that 
coordinate rescaling is an isometry of AdS. QCD is ap- 
proximately conformal at high energies (due to asymp- 
totic freedom), but not at low energies (due to confine- 
ment). This means that an unmodified AdS^ background 
can not capture the essential features of QCD, and the 
AdS^ space must be modified in some way. Generally, the 
holographic QCD models on the market use an asymp- 
totically AdS^ background to reflect the fact that QCD 
is approximately conformal at high energies. The deep 
bulk region is then modified (that is, cut off in some way) 
in order to model confinement. For example, in hard wall 
models [iHl, , [13, [iZi ; confinement is simulated by cut- 
ting off the AdS space at some finite radius by hand. In 
soft-wall models, a dilaton field is turned in the bulk and 
tuned to smoothly cut off the AdS space to produce a 
linear Regge meson mass spectrum [10, fl7Hig |. as op- 
posed to the quadratic spectrum that is generally seen 
in hard-wall models. Other models take into account the 
back-reaction of the bulk fields on the metric, which can 
dynamically cut off the AdS space [3, [l5| . 

A well-known general theoretical issue with all such 
models is that in the regime where QCD is approximately 
conformal, it is weakly coupled due to asymptotic free- 
dom. As a result, it is far from obvious that the use of a 
classical, weakly curved 5D background is justified, since 
one generally expects that a systematic holographic dual 
of QCD should be a string theory on some highly curved 
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space (as noted for instance in ref. (lO'l)- However, despite 
this possible problem, it is interesting to try to investi- 
gate whether a classical 5D background might serve as 
a phenomenologically useful ad hoc approximation to a 
holographic dual of QCD, and many bottom- up models 
appear to show a remarkable agreement with experimen- 
tal data. 

The AdS/CFT dictionary [H [H] is used in bottom-up 
models to dictate the bulk field content: a p-form 4D 
QFT operator O with scaling dimension A corresponds 
to a p-form bulk field with mass = {A — p){A+p — A). 
Conserved currents in the QFT correspond to gauge fields 
in the bulk, with to| = 0. In bottom-up holographic 
models, it is assumed that this dictionary remains valid 
even when one modifies the 5D background away from 
^^55 . It is not obvious that this should be justified, be- 
cause for instance while the scaling dimension of opera- 
tors is well-defined in a CFT, in QCD the scaling dimen- 
sions of most operators receive scale-dependent correc- 
tions proportional to as- That is, most operators in QCD 
have scale-dependent anomalous dimensions. An excep- 
tion is operators that correspond to conserved currents, 
i.e., those associated with vector and axial currents, and 
the stress-energy tensor, whose anomalous dimensions 
are zero. However, in an exploratory spirit holographic 
models of QCD assume that the dictionary remains valid. 

In order to model some observables of interest, one 
must choose a set of 4D operators and corresponding 5D 
fields that probe the relevant physics. Of course, QCD 
has an infinite number of operators with the same quan- 
tum numbers, and due to operator mixing and the lack of 
any obvious suppression scale, an arbitrarily large subset 
of them can contribute to any given process [25|. As an ad 
hoc approximation, to obtain tractable holographic mod- 
els one simply chooses a minimal set of lowest-dimension 
operators to probe the observables one is interested in. 
The parameters of the resulting 5D models are matched 
using the AdS/CFT dictionary to large Nc QCD in the 
UV regime, where the asymptotic freedom of QCD allows 
reliable perturbative calculations. 

Given the many bold assumptions that are necessary to 
construct holographic models of QCD, their phenomeno- 
logical success is remarkable, and may suggest that the 
assumptions are more reliable than might be expected. 
It is natural to wonder if there is anything in these mod- 
els that can test the reliability of the assumptions. In 
this Letter, we show that the ad hoc treatment of scale- 
dependent operators in holographic models of QCD can 
lead to serious phenomenological and theoretical prob- 
lems. For definiteness, we work in the simple hard-wall 
model of ref. 1] for chiral symmetry breaking, and fo- 
cus on the behavior of the quark condensate {qq) and its 
source iriq. We show that to match the large Nc scaling 
of QCD quantities, the operator-field mapping for the 
quark condensate needs to be properly normalized, and 
this normalization affects the value of the quark mass and 
the quark condensate. Next, we argue that the values of 
the quark mass and the chiral condensate obtained in 
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TABLE I: Field content and dictionary of the model. 

the model from matching to data are inconsistent with 
the matching of the 5D models to QCD in the asymp- 
totically UV region, and discuss the implications of the 
analysis for the treatment of scale dependent operators 
in bottom-up holographic models of QCD. 

I. THE MODEL 

Our analysis is in the simple model of ref. [ij, which 
we describe in this section. QCD has an SU{Nf)L x 
SU{Nf)j^ flavor symmetry in the chiral limit, and the 
model focuses on the lowest-dimensional operators im- 
portant to chiral dynamics. These are the left and right 
handed quark currents qLjp.t'^qL and qR'j^t^'qR, where t° 
are the generators of SU{Nf). They are associated with 
the 5D gauge fields Al and using the AdS/CFT dic- 
tionary. We will work in the Nf = 2 limit in this Letter. 

To model chiral symmetry breaking, we include the op- 
erator qRqL, which acquires a vacuum expectation value, 
and associate it with a massive bulk scalar field 4>. (The 
dictionary between qnqL and (j) includes a factor of 1/z 
due to the dimension of qRqL-) The field content of the 
model and the dictionary are given in Table 1. Note that 
in assigning these 5D masses it is assumed that the scal- 
ing dimensions of the operators are purely classical. Note 
that this is not the case for qRqL in QCD, so that this is 
an ad hoc approximation. The model uses an AdS space 
with a hard-wall cutoff in the IR as the holographic 5D 
background. The metric is 

l2 

ds^ = gnydx'^dx'' = -^{-dz^ + dx^dxi), (2) 

where L is the radius of the AdS space, which will be 
set to unity in the rest of this Letter, and z G [0, z„j]. 
Holographic calculations generally use z = e as a UV 
regulator [2^, with e — >■ at the end of the computa- 
tions. We use the convention that Latin indices take 
values 0, 1, 2, 3, and Greek indices take values z, 0, 1, 2, 3. 
This choice of metric corresponds to the assumption that 
QCD, which lives on the z = boundary, remains con- 
formal until confinement suddenly sets in at an energy 
scale of order 1/zm- While this is clearly a drastic over- 
simplification, it does not affect our general conclusions. 

In terms of the bulk fields that we have defined above, 
the classical bulk action can be written as 

j d5xV5Tr||i?</)p+3|0p-^(F2 + F2)|, 

(3) 
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'iAL,,X + iXAR^, Al,r = Al j^t", and 



Ff,^ = d^A^ - d^A^ -i[A^, A^]. The gauge couphng 55 
can be determined from matching the two-point vector 
correlation function to the leading term in the QCD op- 
erator product expansion for the vector two-point func- 
tion. To do this, one can evaluate the 5D action on 
solutions of the equations of motion of the vector field 
iy = Al + Afi), assuming that the 5D side can be 
treated classically, and take two derivatives with respect 
to the boundary sources to obtain the vector two-point 
function. The vector two-point correlation function is 
/ d^xe'P^ {J^[x)J,{Q)) = {q^q, - Q2g^.)ny(g2), and it 
can be shown that near the boundary, which corresponds 
to large scales in the boundary field theory, one obtains [ij 



ny(g') = -^iog(q2) 



(4) 



where q is the 4D momentum. This can be matched to 
the corresponding QCD expression, coming from a quark 
bubble and valid at high g^, where QCD is weakly cou- 
pled due to asymptotic freedom: 



ny(g') = - 
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log(q2) 



(5) 



Matching the holographic model to QCD implies that 
gl = 127rV7V,. 

To incorporate chiral symmetry breaking, it is nec- 
essary to connect the scalar field 4>, which is dual to 
QrQl, to the chiral condensate {(iRqi) = S = crl and 
the quark mass M = rriql, which is the source of qRqL in 
QCD. There is a straightforward prescription for this in 
AdS/CFT. As shown by Klebanov and Wittenf^^, the 
near-boundary solution of the equations of motion for a 
field (/)'-' dual to an operator O with dimension A is given 
by 



<^°(z,<z) = z 



4-A, 



(9) 



2A- 
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where 4i^{q) is the source for O, and (O) = 
5\ogZ/5(j)Q{q) is the one-point correlation function. 

Before applying this to our problem, we must address 
a slightly subtle normalization issue in the association of 
operators with fields in holographic models. In a field 
theory with an operator O and a source J, one always 
has the trivial freedom to redefine O — >■ aO and J — >■ J/a, 
so that JO is unchanged. In holographic models defined 
by actions such as Eq. [1] this amounts to the freedom 
to take O — > aO and (p^ — > (/)^/a[23|- Many treatments 
of holographic models of QCD implicitly take a — 1. 
However, as we are about to show, this is not generally 
correct, and the issue of operator/source normalization 
turns out to have significant implications for holographic 
models such as the one in this Letter. 

We specialize the discussion above to our case with 
^ = and A = 3, keeping the normalization pa- 

rameter a explicit in the equations. Using the dictionary 



lim^_>o f '/'(•z, q) = aM, we can write the q 
solution for 6 as 
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where M = m^l and E is S = {qLqR) = crl- 

It is now clear how the seemingly trivial issue of source 
and operator normalization becomes important. If we 
take a = 1, as was done in ref. [H implicitly, then Eq. [7] 
has incommensurate Nc scaling: the first term scales as 
N^, while the second scales as N^. It is also not hard 
to show that with a = 1 in this model, the p/ai mass 
splitting scales nI^^ , in direct contradiction with large 
Nc QCD. To be consistent with large N^, QCD, we must 

1 /2 

take a ~ Nc rather than a = 1, so that the p/ai mass 
splitting scales as iV°[2l[ . Of course, this should not 
be surprising: in large Nc QCD, the normalization of 
operators explicitly depends on Nc, so the same must be 
true in the holographic approach. 

To complete the description of the holographic model, 
one must specify the values of a, niq and a. The chiral 
condensate a comes from IR physics, and in this model 
it can be taken to be an input parameter. In fact, the 
structure of the model fixes aja and am^. To determine 
cr/a, one can use the relation[l| 



J TT 



1 dA{z)\ 
dz 



9l 



(8) 



where A(z) is the transverse part of the axial vector cur- 
rent (A = Al + Ar) at g = 0, which satisfies the equation 
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[z-^d,A{z)) - ^-^^A{z) = , (9) 

which is shown in the chiral limit niq = for simplic- 
ity. Since this model obeys the Gell-Mann-Oakes-Renner 
(GOR) relation^ mj/^ = 2mqa, it is possible to calcu- 
late amq and a /a by fitting and to data. 

To complete the specification of the model, it is nec- 
essary to determine the normalization parameter a. We 
will compute a by matching a to QCD explicitly, as was 
first done in ref. [31] , and show that this leads to the cor- 
rect Nc scaling. To do this, we will compare the two-point 
functions for qRqL in the holographic models and in the 
asymptotic regime in QCD, using an identical procedure 
to the one that was followed for in the vector sector. 
We will see that this brings up a troubling issue involving 
the matching of the holographic model to QCD. 



II. MATCHING OF THE SCALAR SECTION TO 

QCD 

The equation of motion of the scalar field for general 
4D momentum q is 

d,{z-'d,q^iz, q)) + ^^ (biz, q) = 0, (10) 
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The sohition involves Bessel functions, and can be 
matched for small q to Eq. [T] To compute the two point 
correlation function, we evaluate the scalar field action on 
a solution (pd of Eq. 1101 which leaves a boundary term: 

S[M = J d^x {z-^M^, q)dAcx{z. q)) ■ (11) 

To find the two point correlation function, we can now 
take two functional derivatives with respect to the source 
VmVz^Q ^(p{z,q) — aM, and in the large q limit we find 
that 

J d'xe'P^ mx)qqm - y g2 (^2^2) ^ ^ _ ^ (^2) 

where we have suppressed contact terms. This can be 
compared to the QCD result for large euclidean momen- 
tum q^ and renormalization scale /i: 

/N 
d^xe^^-{qq{x)qqm^^qHog{q^/pi^) + ... . (13) 

This implies that a = y/Wc/i-K, which matches the Nc 
scaling we expected on general grounds. This identifica- 
tion does not depend on any matching between /i and L, 
since both of them can be reabsorbed into the subleading 
terms in the equations above. Following the fitting pro- 
cedure of ref. 1], but using a as extracted above yields 
ruq = 8.3 MeV and a = (213 MeV)^, which differs from 
the results one obtains with a = 1, where = 2. 29 MeV 
and a = (327 MeV)^ !]. 

Of course, even though niq and a are determined by 
matching to experimental data, they have a very spe- 
cific interpretation in the holographic model that can be 
checked against QCD: M = m^l is the source of the 
QCD operator {(iRqi)- 

This presents a serious theoretical problem, since in 
QCD the value of ruq and a are scale dependent. How- 
ever, the quantities aniq and a /a in the holographic 
model are fixed by matching fitting to the GOR relation 
and relating and a. The computation of a above then 
determines rUq and cr, and since a does not depend on 
a renormalization scale, neither do and cr, in confiict 
with the identification of mq as the source of the QCD 
operator qnqh- 

This situation is actually rather common in phe- 
nomenological models of QCD, where one probes various 
observables that are scale-dependent in QCD in models 
where they do not depend on the scale in any systematic 
way. Examples of this include the treatment of structure 
functions in bag models, and the treatment of chiral con- 
densates in Nambu-Jona-Lasinio (NJL) models[30|. The 
phenomenological models are generally taken to be at 
some "natural scale," generally about 1 GeV , and the 
observables computed in the models are taken to corre- 
spond to QCD quantities evaluated at that scale. This is 
the interpretation taken in ref. Q for rUq and a. 

However, while this phenomenological approach may 
be reasonable in models such as NJL model, it is not con- 
sistent with the general structure of holographic models 



of QCD. One of the greatest attractions of holographic 
models of QCD is that they can be matched to QCD. The 
features and parameters of the holographic models are 
generally taken from the AdS/CFT dictionary, so that 
the 5D model is matched to QCD on the AdS bound- 
ary. To match the parameters of the holographic models 
to QCD, it is necessary to compute correlation functions 
in the bulk, and then match them on the AdS bound- 
ary at high 4D momentum q^ to the equivalent QCD 
parameters evaluated at the same high q^ and renormal- 
ization scale /i^. However, in the matching procedure for 
a above, which is used to find mq, the renormalization 
scale /i does not appear explicitly on the 5D side, since 
the 5D model is classical, so that the scale to which the 
holographic computations are supposed to correspond to 
is not obvious. 

Of course, in holographic models it is generally as- 
sumed that 1/z plays the role of /i. However, this does 
not resolve the issues with scale dependent quantities like 
TUq and cr. The issue is due to the fact that the AdS/CFT 
dictionary relates fields on the AdS z = boundary to 
QCD quantities in the UV. The matching to QCD is done 
at asymptotically high scales, where it is weakly cou- 
pled, and /i — oo, in accordance with the identification 
1/z ^ fi on the holographic side. However, in QCD[29|, 
as /X — >■ oo the quark mass rUq runs to zero and chiral con- 
densate a runs to infinity. Consistency with QCD then 
implies that in the holographic model m^, which is fixed 
on the 2 = boundary, should also be zero, while a in 
the model must diverge because of the GOR relation. 

This is clearly inconsistent with the phenomenology of 
the model, which requires that cr ^ cx) in order to have 
a finite splitting between the p and the ai mesons. We 
note, moreover, that this problem does not go away in 
the chiral limit of = 0, since cr still diverges. 



III. DISCUSSION 

Our analysis above should not be surprising: the con- 
struction of the holographic model required a number of 
ad hoc assumptions that are clearly connected to this 
issue. The behavior of scale-dependent quantities like 
the chiral condensate is an explicit probe of the self- 
consistency of the assumptions. Clearly, if one wants 
to match the key features of large Nc QCD in a consis- 
tent way, it is essential to capture the scale dependence 
of QCD on the 5D side of the model. This amounts to 
trying to improve on the the ad hoc approximations in- 
volved in the construction of the 5D model. 

It is well known that the {qRqL) has a scale-dependent 
anomalous dimension (5(/i) proportional to as to lead- 
ing order. Presumably, incorporating such an an anoma- 
lous dimension would change the 5D mass of the scalar 
field. Since the anomalous dimension depends on the 
scale ^ 1/z, it is reasonable to make the 5D mass in 
a holographic model depend on z, m\ = m|(z), with the 
constraint that lim^^o ^ti^z) = —3. Matching the 5D 
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model to QCD for would then amount to demanding 
that nriq in the holographic model obey the same renor- 
malization group equation as the quark mass in 
The anomalous dimension is proportional to the running 
coupling a^, suggesting that one must also modify the 
5D background to allow the as to run[3]- Also, one can 
include a more general potential in the action for the field 
(^^32,], at the price of increasing the number of parameters 
in the model. 

Although we have focused our analysis on the simple 
model of ref. [l| , the problems with scale dependence ap- 
ply rather broadly to bottom-up holographic models of 
QCD, which as yet have not treated scale-dependence of 
quantities like mq and a consistently. It is an open ques- 
tion as to whether position-dependent 5D masses and 
more realistic 5D geometries could make the treatment 
of the chiral condensate in holographic models consistent 
with QCD. 

The analysis here explicitly demonstrates that po- 
tentially serious consequences arise for scale-dependent 
quantities when one attempts to match holographic mod- 
els to QCD in the weakly coupled asymptotic region. One 
can try to avoid the issues of matching to QCD in its 
weakly coupled regime by imposing a UV cutoff on the 
holographic models, as for instance in refs. [H,!!!]. How- 
ever, refs. [H, 113] showed this that a fit to data with the 
UV cutoff as an additional parameter then gives models 
that are defined on a rather short slice of AdS space, rais- 
ing questions as to whether it is still reasonable to use 
the AdS/CFT dictionary directly. 

Since the reasons for worrying about the consistency 
of the matching to the weakly coupled region are rather 



general, it is not implausible that problems may also arise 
even for conserved currents. Specifically, there are some 
subtleties associated with vector currents that we will 
discuss in a forthcoming publication. 

In holographic models that do not systematically deal 
with scale-dependence, the natural way to evade the 
problems discussed in this Letter is to give up on match- 
ing the holographic models to QCD in the asymptotic 
regime, and simply fit all of the parameters of the 5D 
models phenomenologically. In the analysis above, this 
would correspond to fitting aniq and a /ato data, without 
computing a separately, and giving up on the identifica- 
tion of ruq as the source of the QCD operator qRQL- This 
has the disadvantage of losing many of the theoretical 
connections to QCD. This scenario is not ideal, but it is 
not obviously inconsistent. While it is conceivable that 
it may be possible to treat scale-dependence consistently 
by modifying the holographic models, it is clear that the 
issues that we have discussed in this Letter must be ad- 
dressed in bottom- up holographic models of QCD. 
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